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1 Introduction 



Let M be a compact m dimensional Riemannian manifold with smooth 
boundary DM. Let D be an operator of Laplace type on the space of smooth 
sections to a vector bundle V. Let Dg be the realization of D with respect 
o to the boundary conditions denned by a suitable local boundary operator £>; 

we assume D& is self adjoint. In this paper, we shall study the heat trace 
asymptotics. 

We begin by reviewing the situation in the classical setting and refer 
to |b§, pi], |30[ for further details. We suppose given a decomposition of 



dM = Cn^Cd as the disjoint union of two closed (possibly empty) sets. On 
DM, let u- m be the covariant derivative of u with respect to the inward unit 



normal; we use the natural connection defined by D - see []T(J for details. 
Let the boundary operator 

Bu := u\ Cd © (u- m + Su)\ Cn 

define Dirichlet boundary conditions on Co and Robin boundary conditions 
on Cat. Let be the initial temperature distribution. The subsequent heat 
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temperature distribution u := e~* Ds is denned by the equations: 

(d t + D)u = 0, w(0; x) = and Bu = 0. 

Let {(pi, Xi} be a discrete spectral resolution of Dq] the </>i are smooth sections 
of V which form a complete orthonormal basis for L 2 (V) so that 

B(f>i = and D<pi = X^. 

The fundamental solution of the heat equation is trace class. We define 

a(D,B)(t) := Tr L 2e~ tDB = Eie~ tXi . 



Theorem 1 Let D Co = 0- The heat trace a(D,B)(t) has a complete 
asymptotic expansion as 1 1 of the form: 

a(D,B)(t) ~En>o^A£)t (,wn)/2 . 

TTie asymptotic coefficients are locally computable as the integral of smooth 
local invariants: 

a n (D,B) = f M a n (x,D)dx + J CN a+(y,D,B)dy 
+ Ic D o.~{y,D,B)dy. 

These invariants have been computed for n < 5, see for example 0, |5], 19. 
PB] . . There exists a canonical connection V and a canonical endomorphism 
E so that .D = — Tr(V 2 + E). If D is the Laplacian on p forms, then V is 
the Levi-Civita connection and E is given in terms of the Riemann curvature 
tensor R by the Weitzenbock formulas. Let indices i, j range from 1 through 
m and index a local orthonormal frame {ej} for the tangent bundle of M. 
Near the boundary we choose an orthonormal frame so e m is the inward 
unit normal; let indices a, b range from 1 through m — 1 and index the 
induced frame for the tangent bundle of the boundary. Let L^ 1 be the 
second fundamental form of DM C M. We adopt the Einstein convention 
and sum over repeated indices. 

Theorem 2 

1. a (D,B) = (An)- m / 2 f M Tr(I v ). 
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2. ai{D, B) = (4 7 r)( 1 -)/ 2 i{ J Cn Tr (I v ) - J Cjj Tr (I v )}. 

3. a 2 (D, B) = (47r)—/2 1 { j m Tr [R tJ3l l y + 6E ) + J Cn Tr (2L aa I v + 125) 

+ / CD Tr(2L aa / y )}. 

^. a 3 (D, B) = -(4vr)( 1 - m )/ 2 ^{ / Cjv Tr (9QE + 965L aa + 192S 2 
+ (16Rijji — 8R amma + 13L aa L bb + 2L ab L ab )I v ) 
+ Jc* D ^ (96.E + (16i2i^j — SR amma + 7L aa L bb — 10L ab L ab )I v )} . 

In the classical setting, H C^j is empty so the Neumann and Dirich- 
let components do not overlap. There are, however, physically reasonable 
settings where S := (TdHCtv is a non-empty smooth submanifold of dM 
of dimension m — 2. Drop a solid ball at initial temperature into icewa- 
ter. Supposing it floats, the part of the boundary of the ball which is in 
air satisfies Neumann conditions and the part underwater satisfies Dirichlet 
conditions. Here, B is defined by complementary spherical caps about the 
north and south poles of the ball which intersect in a circle of latitude. 

The setting where S is not empty is known in the literature as the N/ D 
problem. It has been investigated extensively from the functional analytic 



point of view |20| , P5| , |26], |29[. However, there are only some preliminary 
results 0, ^, U available concerning the heat trace asymptotics. It is natural 
to conjecture that Theorem [I] can be generalized to this setting by adding 
an extra integral over E of some suitably chosen local invariant. The point 
of this note is to indicate that the situation is not quite so simple. More 
specifically, we will show that the following conjecture is false. 

Conjecture 3 Let Cn^Cd be a smooth hypersurface in dM. The heat trace 
a(D, B)(t) has a complete asymptotic expansion as t j of the form: 

a(D,B)(t) ~En>oan(A£)t("- m )/ 2 . 

The asymptotic coefficients are locally computable as the integral of smooth 
local invariants: 

a n (D,B) = f M a n (x,D)dx + f CN a+(y,D,B)dy 

+ Ic D a n (v, A B) d v + is A B)dz. 
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Here is a brief outline to the paper. We shall suppose that Conjecture ^ 
holds and argue for a contradiction. In §^|, we discuss some of the functorial 
properties which the invariants would have. In §|3|, we use the local index 
formula in a specific situation to show Conjecture |3| is false at the 03 level. In 
§f| we present some results using a perturbation expansion around an exactly 
soluble, but restricted, case which also relate to this question. In §| we 
conclude by suggesting an alternative form that the expansion might take. 



2 Properties of the local invariants 

We use dimensional analysis to study these invariants; this involves studying 
the behavior of the heat trace under rescaling. Let Q be the curvature tensor 
of the connection determined by D. We assign weight 2 to the tensors R, 
Q, and E. We assign weight 1 to the tensors S and L. We increase the 
weight by 1 for each explicit covariant derivative which appears. It then 
follows that the integrands which can be used to compute the invariants 
{aff , a^, a^) are universal polynomials which are weighted homogeneous of 
degrees (n, n — 1, n — 2). For example, Theorem |2| expresses 0,2 in terms of an 
interior integral of Tr (RijjJv + 6E) and boundary integrals of Tr (2L aa ) and 
Tr (2L aa + 12S); these expressions have weights (2, 1, 1). Thus we see that 

Oq =0, a x — 0, and = c dim(V). 

Calculations of Avramidi 0, [| and of Dowker |Bj suggest that 

c = -f(47r)- m / 2 . 

We shall work with the coefficient to show Conjecture |3| fails. Let 
L s be the second fundamental form of X C CV Any invariant which is 
homogeneous of weight 1 can be expressed linearly in terms of the tensors 
{S, L s , L dM }. Thus, in particular, the geometry of the operator D does not 
enter into af. We choose a local frame field on E so that e m -\ is the unit 
normal of £ C Cn and so that e m is the unit normal of Cn C M. Thus the 
structure group is 0(m — 2). Let 1 < u < m — 2. We use H. Weyl's theorem 



31| on the invariants of the orthogonal group to see there exist universal 



constants q so that 

ag(z, D, B) = Tr {( Cl L* u + c 2 L™ + c.L^My + c 4 S}. 
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Let M = Mi xM 2 where M 2 is a closed manifold. Let D := D^l + l®-^ 
where the operators -Dj are operators of Laplace type over Mj. We use a 
suitable boundary condition B\ for .Di to induce a corresponding boundary 
condition B for D. The discrete spectral resolution for (D, B) is given by the 
tensor product of the corresponding discrete spectral resolutions for (D\, B\) 
and D 2 . Therefore 

a(D, B){t) = a(A, B x ){t) ■ a{D 2 ){t) so (2) 

an(y, D ,B) = J2 P + q =nap(yi, Di, B^a^, D 2 ), and 
a^O, -D, £) = E P+q=n dp(zi, Di, B x )a ? (x2, £> 2 )- 

^4 p r i or ^ the constants Cj of equation (|l|) could depend on the dimension 
m but the usual trick of dimension shifting using equation (|2|) and taking 
product with a circle shows the constants q are dimension free modulo a 
multiplicative normalizing factor involving suitable powers of 4n. 

We say that an operator A : C°°(Vi) — ► C°°(y 2 ) is of Dirac type if the 
associated second order operators D 1 := A* A and D 2 := AA* are of Laplace 
type. We assume given boundary conditions B % so that A intertwines the 
spectral resolutions of and Dg 2 - We define: 

index(A) := dim ker(Z)gi) — dim ker(_Dg 2 ). 

The cancellation argument of Bott then shows 

a m (D\ B 1 ) - a m (D 2 , B 2 ) = index(A) and (3) 
a n (D\ B 1 ) - a n (D 2 , B 2 ) = if n ^ m. 

We shall apply this observation in §|3| to the de Rham complex with abso- 
lute or relative boundary conditions. McKean and Singer |2l| used equation 



(0) to prove the Gauss-Bonnet theorem if m — 2; we refer to JXTJ] for a 
discussion of the higher dimensional setting. 



3 An example on the cylinder 

Let M := [0, 1] x S 2 be the cylinder with the standard product metric. We 
consider the de Rham complex and let A e, ° be the bundle of even and odd 
differential forms. Let D e, ° be the associated Laplacians. Let x G [0, 1] be 
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the normal variable and let G S 2 be the angular variable. We then have 
natural decompositions: 

A M = A | © dx A A% D e M = {-dl + A|) © {-d 2 x + A° 5 ) 

Am = A^®&A A e 5! D° M = {-dl + A° 5 ) © {-dl + A e 5 ). (4) 

If (pi G C°°(A,s) are differential forms on M taking values in A$, then absolute 
and relative boundary conditions are defined by the operators: 

i3 a (0i + dx A 2 ) := {d x {(pi) © <M and 
£ r (0i + dx A 2 ) := {9 X (0 2 ) © 0i}. 

We have 

{d + 5) M {<Pi + dx A(j) 2 ) 

= {-d x <j) 2 + {d + + dx A - (d + 5) 5 2 }- 

Suppose that A0 = A0 on M. If B a (j) = on an open subset O C <9M, then 

4{-<9 x 2 + {d + 5) s <f)i}\o 

= (-A s + \){<fa\o} + (d + 8) s {d x <l> 1 \o} = 0fmd 
{dxfa -{d + <5) 5 <Mlo = ~{d + 5s){0 2 |o} = 

so B a {{d + 5)<p} = on O. Similarly if B r <p = on O, then 

{-<9 X 2 + {d + = ( d + $s){<t>i\o} = and 

d x {d x (f)i - {d + 5) S (/)2}\o 

= {A s - X){(f) 1 \o} ~{d + 6) s {d x <h\o} = 

so B r {{d + S)<f)} = on O. Thus {d + <5) preserves the eigenforms of the 
Laplacian with either absolute or relative boundary conditions. 

Let B~ denote pure Dirichlet and B + pure Neumann boundary conditions. 
The structures decouple and we may decompose 

(Dm, B a ) = {-dl + A s , B+) © {-dl + A° s , B~) 
{D M , B a ) = {-dl + A s , Br) © {-dl + A°, B + ) 
{D M , B r ) = {-dl + A s , B~) © {-d 2 x + A s , B + ) 
(£>£,, B r ) = {-dl + A s , B + ) © {-dl + A% B~). 
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The interior invariants vanish if n is odd. Thus for dimensional reasons, 



a 3 (D e M )(x, 0) = and a 3 {D° M )(x, 6) = 0. (5) 
We use Theorem (2) to see that 

af(y,-d 2 x ,B ± ) = ±\ioTyed[0,l}. 

The Dirichlet and Neumann boundary conditions for A decouple on S 2 . Since 
the structures on M are product, we may apply equation (H) with n = 3, 
p — 1, and q = 2 to the decompositions of display to compute: 

af M (Z^, 9) = J{af(9, A|) - af (9, A° s )}, 

a 9 3 M (D° M ,B a )(y,e) = J{-af (9, A|) + af (6, A|)}, 
af M (Z^,23 r )(y,9) = ±{-af (6, A|) + of (9, Ag)}, 
af M (Z^,£ r )(y,9) = ±{af (6, A|) - af (9, A° 5 )} 



for y G <9{[0, 1]}. McKean and Singer [2~I] showed that 



af(6,A|)-af(e,A° 5 ; 



2tt' 



this also follows from Theorem [| since i? = for A| while Tr E = —IR^ji 
for A£. Consequently, we compute the index densities: 

{af\D^ B a ) - a 9 3 M (D° M , B a )}(y, 9) = £ 

{«r B r ) -af M (I^,£ r )}(y,9) = -£ (6) 

Let C a; 5 C 5 2 and CV 5 s C S 2 be complementary spherical caps about the 
north and south pole in S 2 . Let C a> M := {0, l}xC fli s and CV 5 m '■= {0,l}xC nS 
give a corresponding decomposition of the boundary of M. Let B be the 
boundary condition B a on C U) m and £> r on C^a/. The decompositions of 
display ([|) induces corresponding decompositions of B as the sum of two 
boundary conditions of the form we have been considering: 

Wi) = Hc aM © (dx<t>i)\c r ,M and 
i3(dx A 0i) = (d^) | CoiM © 0i|c,, M 

As the metric is product, L®^ 1 = 0. Since S = 0, the only non-zero term in 
given in equation ([!]) is CiL„ M Tr (/). This term is sensitive to the normal of 
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£ C 8M. When studying D e M or D° M we have D/N boundary conditions on 
half the bundle and N/D boundary conditions on the other half the bundle. 
Thus this term cancels and we have 

a*(D e M ,B)(x,e) = and af (D° M , B)(x, 9) = 0. (7) 

We may therefore use equation equation (||), and equation (^) to see 

a 3 (D e M , B) - a 3 {D° M , B) = ^{ V0 l(C o , M ) - vol(C r>w )}. 

Thus this difference is not an integer if the spherical caps are not hemispheres. 
On the other hand, since d + 5 intertwines the eigenvalues of the Laplacian 
with either absolute or relative boundary conditions, we can use equation (||) 
to see that 

a 3 (D e M , B) - a 3 (D° M , B) = index(d + 5,B) G TL. 
This contradiction shows that Conjecture |3| is false. 



4 Special cases and perturbative expansions 

The special case of the N/D wedge has been considered locally by Avramidi 
and globally by Dowker || , where the N/ D hemisphere problem was also 
introduced. In this section we first enlarge on this last example for which 
the N/ D problem can be solved explicitly in terms of known functions, and 
then perturb it to give a more general geometry. 

Consider the hemisphere placed at y > such that its boundary is at 
(p = and (p — 7T with 9 6 [0, 7r]. The metric in the standard (6, (p) polar 
coordinates and the Laplacian are given by: 



ds 2 = d9 2 + sin 2 6dip 2 and A = — — — — — sin^- 



d* 1_ d_ .^ e _d_ 

sin 2 9 d<p 2 sin 9 89 89' 

We are interested in the eigenfunctions Y(ip, 9) which satisfy the N/D bound- 
ary condition. To get an idea of the crucial difference between these condi- 
tions and the classical Dirichlet and Neumann ones, we will provide details 
for all of them. 

The starting point is the separation of variables, Y(<p,6) = $(^>)H(0), 
which leads to the usual differential equations: 

= $"(^)+ y u 2 $(^) 

= sin 9 {sin 9 E' {9))' + {X 2 sin 2 9 - fi 2 )E(9). 
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With the substitution £ = cos 9, E(9) = one finds 

o = (l-a^o-w + ^-^j^o- 

The general solution for <3> is $(y?) = asin(fup) + bcos([Mp), with a, b G (D. The 
equation for u is the differential equation of the associated Legendre func- 
tions. Which eigenf unctions survive as being linearly independent depends 
on the value of fi and thus on the boundary condition. We refer to |27| for 
further details. 

4.1 Dirichlet boundary conditions mean that $(0) = $(7r) = and hence 
ji G IN. Independent solutions are -P^(£) and Qf (£) with I > fM, I G IN, and 
eigenvalues A 2 = Z(7 + l). Imposing square integrability on the eigenfunctions 
shows that the discrete spectral resolution is given by the functions: 

Y{<p,0) = A/lsin(/i^)P; M (cos^), fi G IN, l = n + fi, n G 1N , (8) 

with M\ a normalization constant. For 9 — > and — ► 7r, which is the limit 
to the north and south poles, we see that 

P/\cos0) ~ (1 - cos 2 #f /2 = (sintff 

and so the eigenfunctions are differentiable at the poles, which of course, 
must be the case. We refer to [ 15] for further details. 

4.2 Neumann boundary conditions mean that $'(0) = $'(7r) = 0. This 
yields the quantization condition, /z G IN and, arguing as before, we see that 
the discrete spectral resolution is given by the functions: 

Y((p,9) = Af 2 cos (n^Pf (cos 9), /i G 1N , I = n + /i, n G IN . 

Again, these are differentiable everywhere, including the poles. 

4.3 N/D boundary conditions mean that $(0) = $'(77) = 0. Conse- 
quently = n + 1/2, n G 1N , or, stated differently, \i = TO/2, m = 1,3, 5, ... 
Thus the discrete spectral resolution is given by the functions: 

Y(<p,9) = A/" 3 sin(/^)Pf M (cosfl), fx = f , to = 1,3, 5,..., n G IN , 

with I = n + fi. The vital difference is that now /i is a half-integer for which 
the limiting behaviour of the eigenfunctions near the poles, i.e. the edge, is 

P[ fh/ \cos9) ~ (sin^ 2 , 
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so that the eigenfunctions are not differentiable at the edge. This is the 
crucial difference between the N/ D problem and the classical problems, and 
is responsible for the non-standard small-t behaviour of the heat-trace. 

Having given the eigenfunctions for the hemisphere, for which the bound- 
ary extrinsic curvature vanishes, we now provide a perturbation approach 
which allows account to be taken of the influence of an extrinsic curvature 
at a boundary. Some general developments are given first and applied to the 
hemisphere later. 

Assume the unperturbed situation A0a = A 2 0a with 4>\ dM = with non- 
degenerate eigenvalues. Parametrise the boundary by y and take s(y) to be 
the geodesic distance from it. We define the perturbed boundary dM. e by the 
function es(y) with e very small and call the resulting manifold Ai e - Then, 
to order e, the perturbative formulation of the problem reads as follows, 

Ai/j a = a 2 ip a , ip a \ dMe = °> 

where the initial ansatz for the eigenfunctions and eigenvalues is given by: 

ipa = 0a + e0„, a 2 = \ 2 + er] a . 

The perturbation of the eigenvalues, 7] a , is determined by: 

with the exterior normal d n ; see, for example, [B^]. For Dirichlet conditions 
one can use the identity (p'\dM e — ~\4>\\dM t together with the expansion 

= -es{y)d n (j)x\ dM + ... to see (9) 
rjcc = -~ JdMc 0a<9„0 a = - f dMe fadnfi S dM s(y)\d n (p x \ 2 dy. (10) 

It is important to note that the expansion (||) is well defined for the Dirichlet 
and Neumann eigenfunctions (§). 

In the case of degenerate eigenvalues, which is needed here, the situation 
is slightly more complicated. Let j index the degeneracy. In this case one 
can obtain the secular equation 

X j c id tfJ ki + B%) = 0, with (11) 
B& = - JaMfa<ftHv)M)- (12) 
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Similarly, when considering Neumann conditions and when perturbing 
the shape of the boundary, the equation analogous to (|TTD becomes 



E i c ji (^ fcj - + Sg) = 0, with 
B^^IaM^siy^ji). (13) 

Some of these formulae are implicit in the work of Frohlich, ||, but can 
be traced back to Rayleigh, p8| . The extensive discussion in Morse and 



Feshbach, [p2fl , contains all that one needs. A more mathematical treatment 
is given by Garabedian and Schiffer, || Chap.V, based on Green's theorem 
and Hadamard's formula. 

We apply these developments to the hemisphere. It seems natural to dis- 
place the entire boundary by an amount, say e, perpendicular to the equator, 
i.e. the rim, thus making the new manifold a cap. Geometry shows that the 
extrinsic curvature, L dM , of the perturbed boundary equals e. Then the 
relevant integrals for the eigenvalue perturbations are, 

B D m , = -2Af?^>L 9 M tf^P? ( C os6) Pf' (cos 6), and (14) 
= -2AT 2 V 2 L^ JZ-$- g P?(cos6)P?'(cos6). 

The behaviour of the Legendre functions for 9 — > 0, 7r shows that these 
integrals exist and, although an explicit evaluation is tedious (they are sur- 
prisingly not listed in standard references like Jl], [13|]), in principle this deter- 
mines the eigenvalues to order e and undoubtedly would reproduce the correct 
leading heat kernel expansion. We have not pursued this calculational check. 

It seems reasonable to apply the same perturbative approach to the N/ D 
problem. In this case the N and D parts contribute additively to the secular 
equation, 

E%(tt, + 5g + 5© = o, 

j 



where the definitions (|i2l) and (|13[ ) still hold but involving the eigenfunctions 
of the N/D problem. As one soon realizes, all integrals exist, except the ones 
with yL — jj! — — 1/2. The reason may be found in the use of (|9|) which cannot 
be applied for these modes because it leads to divergences at the edges of the 
manifold. To avoid the use of (|9p we revert to (|T0D and evaluate just these 
awkward modes on the perturbed boundary, dM. t , directly. 
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We return to Dirichlet conditions, where everything is well defined, to 
illustrate the situation. The boundary dM. e can be parametrised by noting 
that along it, y — e — sin^sin^, i.e. 



sin 



-i ( e 



Vsin 9 



x > 0, and <p = n — sin 



-i 



sin 9 



x < 0. 



To leading order in e, the geometrical quantities (normal derivative and vol- 
ume element) on dA4 £ agree with those on dM. and, up to irrelevant correc- 
tions, one finds 



B 



D 



2,/ T dM 



-2N{n'L 



d9 
sin 9 



PI 1 (cos 9) Pi 1 ' (cos 9) x 



sin 



/i sin 



sin 9 



cos 



/i sin 



sin 9 



(15) 



To the relevant order, the integration limits can be set to and ir and com- 
parison of (0) with (|l^) shows that the expansion with respect to e of the 
eigenfunctions evaluated at dAi e leads to the occurrence of a factor e/ sin^ 
to give agreement with the perturbation form. 

Turning to the N/D problem, we are therefore led to consider integrals 
of the type 



which are well defined in the limit e — > for all values of /z, fi' except fi = 
fi' = 1/2. For these "critical subspaces", using the explicit Gegenbauer 



representation [13 



P;%(cos9) = 
the relevant integrals have the form 



s\n(v0) 



jr e d9 sin2[{n + m . 



which behaves like loge = logL as e — * 0. All other eigenfunctions yield 
an 0(e) term to leading order and will not change this log behaviour. These 
remarks suggest the existence of a term L dM log L 9M in the heat trace ex- 
pansion which is an indirect indication, via dimensional arguments, of the 
appearance of logt terms as well (see § 5). 
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Thinking about higher order perturbation theory, it is expected that 
modes which are differentiable only (k — 1) times, fcelN, lead to the occur- 
rence of e fc loge with associated logt terms in the heat trace expansion. 

We stress that, for the perturbed geometry, the extrinsic curvature L dM 
does not vanish at the edge and this is the cause of the trouble. When the 
boundary perturbation, s(y), vanishes at the edges, the perturbation is well 
defined and we expect the standard trace expansions to hold. 



5 Conclusions 

We have shown that if there is an asymptotic expansion of the heat trace for 
the N/D problem, then it is not as simple as in the standard setting. It is 
possible that logt terms enter, the generic behaviour for singular situations. 
There is one setting where such terms are known to arise. If instead of 
local boundary conditions, spectral conditions are imposed, one has a partial 
asymptotic expansion of the form, 

a(D,B)(t) = a n {D,B)&- m ^ 2 + 0(H). 



Again the invariants are locally computable and we refer to |7J, |12| for formu- 
lae if n < 3. However, the complete asymptotic expansion involves non-local 
and log terms HIT , 18 ]. Thus perhaps Conjecture [3] should be replaced by an 



asymptotic expansion of the form 

a(D, B) (t) ~ En> (<*„( A B) log t + a n (D, B) ) t^"' 1 2 



where the leading term a n (D,B) = J E a^(z, D, B)dz is locally computable 
and where the difference 



a n (D, B) - J M a n (x, D)dx - J Co a n (y, D, B)dy - J Cn a+(y, D, B)dy 

is a non-local invariant determined by the behavior of D and B near S. As we 
have argued in Section 4, logt terms may occur to compensate the \ogL dM 
terms present in the perturbed hemisphere example. Further study of the 
heat trace asymptotics of the D/N problem seems indicated. 
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